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Abstract 



We develop a novel eikonal expansion for the Cooperon to study the effect 
of space- and time-dependent electric fields on the dephasing rate of disor- 
dered conductors. For randomly fluctuating fields with arbitrary covariance 
we derive a general expression for the dephasing rate which is free of infrared 
divergencies in reduced dimensions. For time-dependent external fields with 
finite wavelength and sufficiently small amplitude we show that the dephasing 
rate is proportional to the square root of the electromagnetic power coupled 
into the system, in agreement with data by Wang and Lindelof [Phys. Rev. 

Lett. 59, 1156 (1987)]. 
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I. INTRODUCTION 

The dephasing time r^ of a particle in a quantum system is the time over which its wave- 
function maintains phase coherence. For diffusive systems in dimensions d = 1 and d = 2 
the dephasing time can be obtained experimentally from the weak-localization correction to 
the conductivity in a magnetic field. Recently Mohanty, Jariwala and Webb |IJ used this 
method a measure r^ as a function of the temperature T in gold wires. They found that at 
sufficiently low temperatures r^ approaches a finite value. This is in disagreement with the 
generally accepted point of view H] that in thermal equilibrium r^ should diverge as T~ p , 
p > 0, for T — > 0. A proposal that the observed saturation of r^ is an intrinsic ground-state 
property of interacting electrons in a random potential |§ has been heavily criticized [|J. 
Until now, there is no general agreement on the microscopic mechanism for the observed 
saturation of r^ at low temperatures. 

One possible reason for this saturation might be external microwave radiation that is 
unintentionally coupled into the system [Q. According to Ref. || the dephasing rate 1/tac 
due to a time-dependent but spatially constant electric field is proportional to the microwave 
power P absorbed by the system for small P, and crosses over to a P 1 ' 5 -law for larger P. 
Experimentally the effect of microwaves on weak localization has been studied by several 
authors PHID!- Most of the data by Wang and Lindelof 0] can be fitted with 1/tac cx P 1 ^ 2 - 
Below we shall offer a simple explanation for this behavior. More recent data by Webb et al. 
T0|1 are consistent with a P 1//5 -law at high powers P, but in this experiment the microwaves 



seem to heat the sample, an effect which has not been taken into account in Ref. ||, and 
which we will neglect as well. 

Following Altshuler et al. @,|ll[] , we define the dephasing rate in the diffusive regime via 



the weak localization correction 5a to the static conductivity. In the presence of an external 
electric field 5a can be written as |]6| JTT|J12] 

a r°° I r T ° 
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ni>d 



f°° 1 f T t t 

/ dt J im T^F dt C(r,r,-,--,t ) , (1) 



where ctq is the Drude conductivity, v& is the d-dimensional density of states, and r e i is the 



momentum relaxation time. The Cooperon C satisfies [12-tt 



dt + DP 2 r + r + i[V{r, t + t)- V{r, t - t)] 



xC(r, r', t, t', t ) = S(t - r')5(t - t') 



(2) 



where P r = — zV r is the momentum operator, D is the diffusion coefficient, and the phe- 
nomenological cutoff T describes dephasing due to processes that are not explicitly treated 
in our calculation, such as inelastic electron-phonon scattering or electron-electron scattering 
with energy-transfers \u\ £, T. The potential V(r,t) is related to the longitudinal electric 
field via eE(r,t) = V T V(r, t), where — e is the charge of the electron. For Tq = V = 
the ^-integration in Eq.(|l]) diverges in d < 2 at the upper limit, but for finite T or V the 
integration is cut off at some finite time, the so-called dephasing time. 

II. EIKONAL EXPANSION 

Formally Eq. (0) looks like the differential equation for the imaginary-time single-particle 
Green's function of an electron in a fluctuating external potential 



V to (r,t)=i[V(r,to + t)-V(r,t -t)]. 



(3) 



The real time version of this problem has been discussed extensively in the quantum field 
theory literature [ 15| . By means of a simple modification of the method developed by E. S. 
Fradkin |TJ] (see also Ref. |Tj|) we obtain the solution of Eq.(0) in the following form fl2 

dk 



C{r,r',t,t , ,t ) = e(t-t') 
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-(r +Dk 2 )(t-t')-F(k,t-t';r,t) 



(4) 



where the function F satisfies the eikonal equation 



d T + dt + D(P 2 r + 2k ■ P r ) F(k, r; r, t) 



V t0 (r,t)+D P r F(k,r;r,t) 



(5) 



with boundary condition F(k, 0; r, t) = 0. Although in general Eq.(|5|) cannot be solved 
exactly, we can easily obtain an expansion of F in powers of the potential V to (r,t). Setting 



F(k,r;r,t) = £F n (k,r;r,t) J 



(6) 



n=l 



where F n involves by definition n powers of Vt , successive terms may be calculated re- 



cursively [0,|T3|,|IB| . For our purpose we need only the first two terms. The linear term 



is 
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Here 



(7) 



V to (q,u) = i[e-^V(q,u) - e**»V(q, -u)} 



(8) 



is the Fourier transform of Vt (r,t), where 



V(q,u) = fdr j dte- i{ * r - ut) V(v,t) 



From Eq.(j|) we find that the quadratic term is related to the linear one via 



F 2 (k,r;r,t) 



The explicit calculation yields 



^■^imi 
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2Dq! • q 2 D(qf + 2k • q 1 + 2qx • q 2 ) - ia>i 

e [D(q 2 .+2k-q2+2q 2 -qi)-iw 2 ]T _ ]_ ~| 
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(10) 



(11) 



Because we have made the diffusion approximation, the momentum-integrations in 
Eqs.(|7|,[n|) are restricted to |q|, |qi|, |q 2 | < l/£ = l/(vpT e \), and the frequency-integrals 
to \u\, |co>i|, |cc^2 1 < 1/ T ei- Here I is the elastic mean free path and Vp is the Fermi velocity. 
For brevity we have not explicitly written out these cutoffs in the above expressions. 

III. DEPHASING DUE TO RANDOM FIELDS 

A. General case 

We now assume that the potential V(r, i) is a random function with zero average and 
general covariance 

(V(q, u)V(q', J)) = (27t) d+1 5( q + q[)6(u + iv')g(q, u) , (12) 

where (...) denotes averaging over the probability distribution of V. The dephasing rate is 
then defined in terms of the average (5a) of Eq.([|). Using the fact that after the averaging 
the Cooperon is independent of the time to, we obtain from Eqs.( 



{Sa ) = -^- [°°dt ( Jk -(ro+^)*-r(M) (13) 

^ 7 7TU d Jr cl J (2%) d ' l ^ 

where 

r(k,t) = -ln(e- F(k ' t;r ' t/2) ). (14) 

We now perform a linked cluster expansion of T(k, t) in powers of the correlator g>(q, uj). To 
first order we find 

r(k,t)^r 1 (k,t) + r 2 (k,t), (15) 

where 

r 1 (k,t) = -^(F 1 2 (k ) t;r ) t/2)), (16) 

r 2 (M) = (F 2 (M;r,t/2)>. (17) 



Introducing the notation -Ek(q) = D{q 2 + 2k • q) we obtain 

dqdio 
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B. Nyquist noise 

As a special case let us assume that the potential V is generated by equilibrium fluctua- 
tions of the electric field due to the thermal motion of the electrons (Nyquist noise). Then 
g(q,u) is determined by the fluctuation-dissipation theorem, which implies [|ll[ 



0(q,w) = -/ q coth(— )Ime 1 (q,u) 



(20) 



where / q is the Fourier transform of the bare Coulomb interaction, and the dielectric function 
of the system is in the diffusive regime given by 

Dq 2 



e(q,w) = l + /q^d 



Dq 2 



ILO 



(21) 



The corresponding dephasing rate in d — 1, 2, 3 was calculated by Altshuler et al. [[TlJ] with 
the help of a Feynman path-integral representation of the solution of Eq.(|2|). It is instructive 



to reproduce the results of Ref. jnj within our eikonal expansion, treating d as a continuous 
parameter. It is important to note that in the derivation of Eq. (Q) the electromagnetic field 
is treated classically |Tl|JT^] , so that our calculation takes only into account low- frequency 



Fourier components of the screened Coulomb interaction, with \u>\ £T. The correlator 
g(q,ui) in Eqs. ([T8| , [T9D can then be approximated by 



IT 
g( q ,uj)^&(T-\uj\)—^, (22) 

Because the momentum integrals in the following analysis will be infrared and ultraviolet 
convergent, no further cutoffs are needed in our approach. 

The dephasing rate is determined by the long-time behavior of T(t) = T(k = 0,t). In 
d < 2 we find for t > 1/T 

Tt 2 ~i 
r(t) ~ C d — -r , (23) 

u d D2 

where the numerical constant C d is given by 

o3— d 
Cd = <,o AM A * ■ (24) 

71-2(2 -d) (4 -d) 
Because T(t) grows for large t faster than linear, the term T Q t in Eq.(]I3]) is negligible, and 
the dephasing rate may be defined by T(r^) = 1, which yields 

C d T 



1 

7> 



u d D d / 2 



(25) 



In d = 2 we find for t > 1/T 



^™. 



so that r^ satisfies 



1 _ ln(Tr )T 



T0 2-ku 2 D 

Keeping in mind that in a good metal 2ttv 2 D 3> 1, this implies to leading order 

1 ln(27rz/ 2J D)T 



(27) 



(28) 



70 2nv 2 D 

in agreement with Ref. ||11| . Finally, in d > 2 the long-time behavior of T(t) is dominated by 
the term proportional to t in the second line of Eq.(|TJ5p. The total dephasing rate in d > 2 
can then be written as 



If dqdu 



f dqdu 



Tj, ' J (27T 



Dq 2 



vr (Dq 2 ) 2 + uo 2 



T + C d T*/(v d Di), (29) 



where 



2 2 ~ d 

Cd = — d : 3— • (30) 

7rl(d-2)r(|)sin(j7r) 

Note that for d — > 2 both prefactors Cd and Cd diverge as 7T -1 |d — 2| -1 , signalling logarithmic 
corrections in d = 2. The term in the braces of Eq.(^) is the dynamic structure factor of 
the diffusing electrons in the regime \u\ £, T, where the detailed balance factor e~ w l T can 
be replaced by unity This term agrees with the semiclassical dephasing rate derived by 



Chakravarty and Schmid |17]], see also Ref fl8" |. It should be kept in mind, however, that 



Eq.(p9"D is only valid in d > 2, where the integral is finite. Diagrammatically Eq.(^) ignores 
vertex corrections in the Bethe-Salpeter equation for the Cooperon, which become important 



in d < 2. In this case Eq.(^) should be replaced by the more general expressions ( |T8| , |T9"| ) , 
which are free of infrared divergencies. 

Our T(t) corresponds precisely to the function fd{t) introduced in a recent preprint by 
Golubev and Zaikin [|TUJ. To make contact with this work, let us assume that Eq.(^) remains 
valid for potentials V(r, t) due to Nyquist noise with frequencies in the range |o;| < l/r e i, 
so that in Eqs. ([l8| , [T9D we may approximate g(q,u) ~ ©(Td 1 _ |^|)^ c oth(^)/(z/ rf Z)q 2 ). 
Then we obtain in the quantum regime Tt ^C 1 in quasi Id to leading order T(t) ~ 
v^/r^V^ t/{v\D^), which agrees with the leading term in the expansion of the func- 



tion fi{t) given in Eq.(28) of Ref. [|T^]. Moreover, also the subleading corrections to /i(t) 



given by Golubev and Zaikin [[HJ can be obtained within our eikonal expansion JT2] . Hence, 
at least in d — 1 we can completely reproduce the behavior of the function fd(t) discussed in 
Ref. |19|] if we assume that the differential equation (|]j for the Cooperon remains valid when 
the classical potential is replaced by a quantum field mediating the Coulomb interaction. This 
replacement has been claimed by Golubev and Zaikin to be consistent with a fully quantum 
mechanical calculation. It is the origin for the discrepancies between their work and Ref. 

8 



[111. 



IV. DEPHASING DUE TO EXTERNAL FIELDS 

In the experiments [0,|1O[ microwaves are coupled into the system via an antenna at- 
tached to a suitable waveguide, such that alternating longitudinal currents are induced in 
the sample. Keeping in mind that the precise way in which the microwaves couple into the 
system is not known, we assume for simplicity that the electrons feel a longitudinal electric 
field of the form 

E(r, t) = E cos(q ■ r - u Q t) , (31) 

where u is the microwave frequency and the wave-vector q depends on the geometry of 
the waveguide and the antenna. The corresponding potential in Eq.(pJ) is 

V(r,t) = V sm(q -r-uj t) , (32) 

with Eq = qoVo/e. Let us emphasize that this is the total screened potential, which is the 
sum of the external potential and the induced potential. Of course, the Maxwell equation 
inside a metal contains a dissipative term, so that the field inside the metal is not given 
by a simple propagating wave [pOR . In general we expect that the field distribution inside 
the metal depends on the boundary conditions and on the precise manner in which the 
microwaves are coupled into the system. Such a calculation is beyond the scope of this 
work. However, if the spatial variation of the field is sufficiently slow (i.e. go is sufficiently 
small) the field inside the metal can still be approximated by a plane wave. To estimate 
the upper limit for go where this approximation is correct, let us assume that the external 
potential applied to the electrons is V ext (r, t) = Vq x1 sin(q ■ r — oj^t). The total potential is 
then 



y(r,i) = y o ext Ime- 1 (q ,w )e ,(qoMot) , (33) 



where the longitudinal electric function in the diffusive regime is given in Eq. (pip . Screening 
can be ignored in the regime where e(qo,u;o) can be approximated by unity. From Eq.(pT|) 
it is easy to see that in d = 1 this is the case when 

go « /f , (34) 

while in d = 2 the external field is effectively not screened if 

* <<: ^ko ■ (35) 

Keeping in mind that in a good metal 27re 2 z/ 2 > l/£, we see that Eqs. ( |34"1 , |35D together with 
u> < 1/t c i are more restrictive than the condition q < l/£ which has to be satisfied in 
order to use the semiclassical Eq.(Q). We assume that q is sufficiently small so that the 



inequalities (|34] , |35D are valid. From Eq.(H) we then obtain 



T/ext 

Vtofau) = ^(2vr) d+1 [5{u - cu ) - 6{u + co )} 



x 



-iuigto 



<5(q-qo) + ^ oto 5(q + q )l • (36) 



This expression is now inserted into the general results given in Eqs. (0) and (|TT|) . For F\ we 
get 

F 1 QL,<r,r,t) = -g L^-^im ^ e ^—e^A 

q { \D(q^ + 2k-q )-tuJo J 

/l _ „-D(q2-2k-q )r p ia;oT \ ~) 

+ e - J (q r- Wo to) Im _^ e ~i^t _ , 37 N 

\D(cti-2k-q )-iLJo J) 

Performing the analogous calculation for the function F 2 given in Eq.([TT|) would result in 
a rather lengthy expression. Since in the following we are only interested in the leading 
terms of an expansion in powers of go > if is more convenient to first expand F\ and then use 
Eq.(|T0p to obtain F 2 . 

If we take the limit qo — ► keeping E constant we obtain a spatially constant field, 
which has been considered in Ref. || . In this limit the microwaves do not affect the density 
of the electrons. For finite qo, however, the electric field induces a density modulation with 
amplitude po — -EWo/(47r). We now show that for sufficiently small microwave power P this 
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leads to a dephasing rate 1/tac proportional to P 1 / 2 , as observed in Ref. f7j. Note that the 
microwave power P coupled into the sample can be estimated as |5|.|10|] 



(EoL) 2 
2-Rtot 



P = ^r^ oc El , (38) 



where L is the effective sample length including leads, and R to t is the total resistance. 

According to Eqs.([I],[|) the weak localization correction 5a is determined by P(k, t; r, t/2). 
In the experiments 0,[nj the wavelength of the microwaves is larger than the size of the 
sample, so that we may expand P(k, t; r, t/2) in powers of go- From Eq.([37]) we obtain 

Fi(k, t;r,t/2) = 2DeE u 2 g 1 (tj t){2k\\sm(tJot ) 

+q cos(cu to)[i(l - 2Dk\t)+2r\\k^} + 0{E Q ql) . (39) 

where k\\ = k • qo/<?o, r\\ = T ■ qo/«?o and 

g^x) = x cos (-) - 2 sin (-) . (40) 



Inserting this result into Eq.fllOD yields 

F 2 (k, t; r, t/2) = AD(eE ) 2 LO 3 g 2 (uj t) sm 2 (oo t ) 

+ O(E 2 q ) , (41) 

with 

g 2 (x) = x + — cos(x) - - sin(x) . (42) 

Note that the second line in Eq.(|39"D is the leading correction for finite go an d small E . 



It can be shown [|12|] that in the limit go — * all higher order terms F n , n > 3 in our 
eikonal expansion vanish, so that the exact solution of the eikonal equation (^) is given 
by F qo=0 = [Fi + F 2 ] qo =o- The fact that the Cooperon in a spatially constant field can be 
calculated exactly is obvious in a gauge where the electric field is represented in terms of 
a vector potential [Bf. It is reassuring to see that our eikonal expansion reproduces the 



exact solution in a different gauge. Substituting Eqs. ( |3"9"| , |4"TD into Eq.(f|) and performing the 
k-integration, we obtain for the weak localization correction to the static conductivity 

11 



6a = a ° <d ' r _^ e -a 9S ( X) - lx 



nu d (4vr J D) o! / 2 Ju T el x d / 2 

I pit e ag 3 (x)cos(2ip)-i/3g 1 (x)cosip 
X - / dtp . . (43) 

nJ o yfl - 2i/3gi(x) cos ip 



where 



2 ' X 



g 3 (x) = x + sin x sin I — ) , (44) 

and we have introduced the dimensionless parameters 

~-^) a (45) 

(46) 
(47) 

Note that after the (^-integration the integral in Eq.fl4"3]) is real, so that we may replace the 
integrand by its real part, which is independent of the sign of E . This is the reason why in 





ul 


p = 


2£|e£ |go 


ul 


7 = 


LU ' 



the definition fl46[) of (3 only the absolute value of Eq appears. Note that (3 is proportional 
to the absolute value of the amplitude p = E q /(4:7i-) of the density wave associated with 
the longitudinal electric field. 

The order of magnitude of the non-equilibrium dephasing rate 1/tac can be estimated 
as 1/tac ~ ^o/xc, where x c is the effective large-x cutoff for the integration in Eq . (|43|) . For 
simplicity we now set 7 = 0. Depending on the ratio a/ (3 = \eE Q \/(2uoqo) = |Vb|/(2u>o)> we 
obtain different behaviors for 1/tac- m the limit a/ (3 ->oo we recover the results of Ref. 
0. In this case 1/tac ^ ^o a f° r a < 1, and 1/tac ^ ojQa l / b for a ~^> 1. Note that a is 
proportional the microwave power P oc _Eq absorbed by the system, while (3 is proportional 
to P 1 / 2 . Obviously, for sufficiently small P we always have /3>«. Then the cutoff for the 
^-integration in Eq.(|4*3]) is determined by (3, and we can set a = to determine x c . The 
^-integration is effectively cut of where the integrand starts to oscillate. Thus we estimate 
x c from the condition (3gi(x c ) ~ 1. For (3 <C 1 the effective cutoff x c is large compared with 
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unity, so that with Eq.(f44|) we obtain l/x c oc (3. This implies 

oc luqP , for a < (3 < 1 . (48) 

tag 

On the other hand, for /3 ^> 1 the x-integration is cut off at small x and we approximate 
gi(x) ~ — Yg. This yields 

f 



oc cj /9 1/3 , for a < /3 and /? > 1 . (49) 

tac 

Hence, for sufficiently small P the dephasing rate 1/tac is proportional to P 1 ' 2 . Moreover, in 
this regime uq/t^q should be independent of ujq as long as the dispersion of the longitudinal 
density wave (i.e. the dependence of go o n ujq) can be neglected. Note that the condition 
P 3> a ^> 7 where 1/tac should exhibit a P 1 / 2 -dependence can also be written as 

Dq 2 » ^r » T . (50) 

ro| 

Thus, for finite T the value of go must be sufficiently large to observe the P 1//2 -law. Keeping 
in mind that according to the conventional point of view @-§[ the intrinsic dephasing rate 
T should vanish for T — > 0, we conclude that at sufficiently low temperatures and small V 
Eq.(p)0"D can be satisfied for experimentally relevant wave- vectors go- Note that according to 
Ref. the observed saturation of the dephasing time, which is typically of the order of a 
few nanoseconds, is due to some external noise. If this is correct, then the experimentally 
observed saturation value of the dephasing rate should not be identified with r . The 
same arguments apply for the condition |Vo| 3> T that follows from Eq.([5(]) together with 
Eqs. ([34]j35|) . For low enough temperatures this inequality should always be satisfied. (Most 



recent experimental data |22|] indicates r — > for T — >• in narrow Ag wires. These probes 
therefore seem to be suiting to test our findings.) In conclusion, in the limit T — > 0, the 
only relevant restriction following from Eq.(|50"D should be \Vq\ <C ujq . 

For convenience, let us summarize the assumptions made in our derivation of the dephas- 
ing time due to external microwave radiation in this section. First of all, we have assumed 
that the microwave field inside the metal can be approximated by a propagating wave with 
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wave-vector q , frequency ujq and amplitude E = q Vo/e. We have argued that this ap- 
proximation is valid in a regime where the dielectric function of the bulk system can be 



approximated by unity, which is possible if go "C Jujq/D in d — 1, and go "C cu /(2ire 2 V2D) 
in d = 2. Throughout this work we have assumed diffusive dynamics and calculated the weak 
localization correction to the conductivity from the semiclassical equation (|2]), which is of 
course only correct for sufficiently small wave- vectors (g < l/£) and frequencies (u> < l/r c {). 
Note that the latter inequality sets the upper limit for the microwave frequency where our 
approach remains valid. If these restrictions are satisfied, then the dephasing rate due to 
microwaves is 

^/Spl/5 for |y Q | > Uq and D ^ Eq) 2 > ^2 ^ 



1 

— oc < 



oj q 2 P for | Vo| > w and D(eE ) 2 < u 2 



tac 



o ■ 

2 



(51) 



(q uJo) 1/3 P 1/6 for \V \ < oj and P>|eP |go > ^ 
qoUJ^VP for |Vo| <^C u and D\eE \q <^ lu 2 , 

where P ex E 2 is the microwave-power absorbed by the system, and Vq = eEo/qo. 

V. COMPARISON WITH EXPERIMENTS AND SUMMARY 

Some time ago Wang and Lindelof have measured 1/tac as a function of P in mag- 
nesium films. Their data from Ref.[7(a)] are reproduced in Fig.|l]. At uj /2tt = 0.66GHz our 
prediction ig> /tac oc P 1 / 2 is in good agreement with the experiment. Although the data at 
uj /2tt = 3.61GHz cannot be fitted by a straight line through the origin, the data roughly 
scale as 1/tac oc Uq 1 for fixed and small P. One should keep in mind, however, that in 
the experiments the precise value of the microwave power coupled into the system was 
not measured, and the power axis for the two sets of data was rescaled differently. Further 
evidence for the P 1 ' 2 -law can be found in Fig. 19 of Ref. [7(b)]. 

Recent measurements of 1/tac by Webb et al. [|10| suggest a P x ' 5 -law for large P in a 
limited range of frequencies. Our calculation shows that the P 1//5 -law should hold as long 
asa> max{/3, 1}, while for j3 3> max{a, 1} we predict 1/tac oc (3 1 / 3 oc P 1 / 6 . Keeping in 
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mind that a/ (3 oc u x we predict a crossover from a P 1 / 5 - via a P 1 / 6 - to a P 1 / 2 -behavior 



as the frequency is increased. The data shown in Fig. 6 of Ref. [1(J are consistent with the 
existence of such a crossover. 

In summary, by means of an eikonal expansion for the Cooperon in a slowly varying 
scalar potential we have derived a general expression for the dephasing rate due to fluctuat- 
ing electric fields in disordered metals. Our method is a physically transparent alternative 



to the path-integral approach used in Ref. [11]. For randomly fluctuating fields with zero av- 



erage and arbitrary covariance the dephasing rate can be obtained from Eqs. ( ]T8|JT9| ) . These 
expressions take vertex corrections into account and remain finite in reduced dimensions, 
where the well-known 0,0 semiclassical result Eq. (|29|) is infrared divergent for some phys- 



ically relevant g(q,u). For example, in the case of 1/ /-noise fllql , where g(q,u) oc 1/u, the 
integral in Eq. (|29D is divergent in d < 2, while our more general result ( [L8l , |l9f) is finite. We 
have also studied dephasing due to external electric fields, and have proposed an explanation 
for the data of Ref. |?j]. Finally, we point out that our eikonal method can also be used to 
calculate the Diffuson in a fluctuating electric field, which is of current interest in several 



contexts 21 
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FIG. 1. Data for loq/tac as function of P 1 ' 2 from Fig. 3 of Ref.[7(a)]. Here P is the microwave 
power coupled into the system. The solid line is a fit of the data at loq/2tt = 0.66GHz to our 
prediction wq/tac °c P 1 ' 2 , which is valid for small P. 



